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The correct expression for the affine parameter 7 along u — uq is given by j{v) — vM[l/{2a), 1 + l/{2a), aA'^v'^" /2], 
and not as appears in the original Letter. A number of quantities which depend on 7 change as a result of this 
correction: For a — 0.3 the correct value for /? is /? = —8.00 ± 0.01, the function f3{a) for supercritical evolutions is 
given by /3(a) oc q,^1-3i±0.07^ ^^^^ ^^^.^ _ _26.9 ± 3.0. The corrected expression for j{v) also changes Figs. 3,4, and 
5. The corrected figures are given below. 
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FIG. 3: Affine parameter at uo, 7,, at which a spacelike singularity is intersected along an outgoing null ray as a function 
of A. Here, a — 0.3. Numerical data are depicted by circles, and the solid line is proportional to A"*"". 



FIG. 4: Affine parameter at uo, 7,, at which a spacelike singularity is intersected along an outgoing null ray u, as a function 
of A. Here, a — 0.1. Numerical data are depicted by circles. 



As the new function 7 is a monotonically increasing function of v, and as 7 — *■ 00 as w — > 00, it satisfies all the 
necessary requirement for comparison between advanced times of different spacetimes. 

The main conclusions of the original paper remain unchanged. In particular, we find that there are initial data 
(subcritical data) for which the singularity has a hybrid structure (consisting of a spacelike sector and a null sector) , 
and initial data (supercritical data) for which the singularity is spacelike throughout. 
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FIG. 5: The exponent (3 as a function of a. The circles are the data points, and the solid curve is proportional to a~^'^^ . 

I thank Asccm Paranjapc for pointing out that there was an error in the function ^{v) in the original Letter. This 
research was supportred by NSF grant No. PHY-9734871. 
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The singularitiy inside a spherical charged black hole, coupled to a spherical, massless scalar 
field is studied numerically. The profile of the characteristic scalar field was taken to be a power 
of advanced time with an exponent a > 0. A critical exponent Ocrit exists. For exponents below 
the critical one (a < ctcrit) the singularity is a union of spacelike and null sectors, as is also the 
case for data with compact support. For exponents greater than the critical one (a > Ocrit) an all- 
encompassing, spacelike singularity evolves, which completely blocks the "tunnel" inside the black 
hole, preventing the use of the black hole as a portal for hyperspace travel. 

PACS numbers: 04.70.Bw, 04.20.Dw 



One of the most intriguing questions in general relativ- 
ity is the inevitability of evolution of singularities under 
very plausible conditions. Indeed, some view this pre- 
diction as evidence that the theory "carries within itself 
the seeds of its own destruction," and look for alter- 
native theories which are singularity free. An alterna- 
tive, more tolerant approach is to view singularities as a 
phenomenon "from which we can derive much valuable 
understanding" Here, we adopt the latter viewpoint. 

There are two ways in which generic singularities in 
general relativity are manifested. The first, and the more 
familiar way, is through unbounded tidal disruptions, like 
the singularity inside a Schwarzschild black hole. Sin- 
gularities of this type are spacelike and deformationally 
strong. Extended physical objects which approach such 
a singularity are unboundedly stretched in one direction 
and compressed in two other directions, such that they 
are inevitably crushed to zero volume. The second way in 
which singularities in general relativity are manifested is 
through breakdown of predictability, as with the Cauchy 
horizon (CH) singularity inside perturbed spinning black 
holes. The CH singularity is non-central, null, and de- 
formationally weak. That is, although it is a genuine 
curvature singularity with infinite tidal forces, the lat- 
ter are (twice) integrable, such that the tidal distortions 
are bounded 0]. This singularity signals an inevitable 
loss of predictability: to the future of the CH singularity 
(if the spacetime manifold can be extended beyond it) 
lurks a region of spacetime which is beyond the domain 
of dependence of any Cauchy hypersurface in the exter- 
nal universe. The occurrence of the null singularity inside 
realistic black holes is of great interest: It leaves open the 
possibility that extended objects could traverse the CH 
only mildly affected, and reemerge in another universe 
(or a distant portion of our universe), in practice using 
the black hole as a portal for hyperspace travel. 

We consider here the characteristic initial value prob- 
lem, in which a pre-existing Reissner-Nordstrom (RN) 
black hole is perturbed nonlincarly by a self-gravitating 
scalar field. The scalar field is specified along a union of 
an outgoing and an ingoing null ray, which is our char- 



acteristic hypersurface. 

The picture of the singularity inside realistic black 
holes which has emerged from analytical and numerical 
studies which assume characteristic data with compact 
support is the following: The onset of the null singularity 
at the CH can be well described by perturbation theory. 
The null generators of the CH focus monotonically, such 
that the singular sphere at the CH shrinks as it moves 
with the speed of light. Eventually the singular sphere 
contracts to zero volume, at which point the causal struc- 
ture of the singularity changes, and it becomes spacelike, 
and deformationally strong. Interestingly, both types of 
singularity coexist inside spherical charged black holes, 
and possibly also inside spinning black holes. This situ- 
ation is described by the lower diagram of Fig. ^ The 
spacelike portion of the singularity inside the black hole 
is not causally connected with the null portion of the sin- 
gularity (except at one point). We can strengthen (in a 
sense to be made precise below) the characteristic data of 
the (nonlinear) perturbation, and thus make the spacelike 
singularity appear sooner. This fits well with the general 
picture: with stronger perturbations the focusing of the 
null generators of the CH is faster, such that the volume 
inside the contracting singular sphere vanishes sooner. 

Can we specify characteristic data which are so strong 
that a spacelike singularity evolves before the onset of the 
null singularity? Below, we shall answer this question in 
the affirmative for a simplified toy model: Perturbations 
slightly weaker than a certain critical value have a sin- 
gularity described by the above picture, i.e., a union of 
spacelike and null sectors. Stronger perturbations involve 
only a spacelike singularity, which completely replaces 
the null singularity and seals off the "tunnel" inside the 
black hole, preventing its use as a portal. 

In this paper we study this question in the context of a 
spherically symmetric black hole, which is endowed with 
electric charge. This is a simplified toy model for a spin- 
ning black hole, as both exhibit similar causal structures 
and instabilities of the inner horizons. This black hole is 
perturbed by a self-gravitating, spherical, massless scalar 
field. Indeed, in this model the above picture was con- 
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FIG. 1: The Penrose diagrams for the two types of singulari- 
ties. Upper diagram: an all-encompassing spacelike singular- 
ity. Lower diagram; The singularity includes both a spacelike 
sector and a null sector at the CH. 

firmed using both numerical simulations and analytical 
analyses for perturbations with compact support ^ Q 
(and one family of characteristic data with noncompact 
support (Sj). (The occurrence of a spacelike singularity 
inside spinning black holes is yet to be confirmed.) Any 
characteristic data with compact support lead to fields 
which decay at late times as an inverse power of advanced 
time along the event horizon, and lead to the null CH sin- 
gularity at early values of retarded times (connected to a 
spacelike singularity at late retarded times). Such data 
are generic, in the sense that they accompany any realis- 
tic gravitational collapse process. The perturbations due 
to the collapse can be thought of as a lower bound on 
the perturbation field. It is of interest to ask what the 
singularity structure is when initial data stronger than 
the lower bound are present. A simple class of noncom- 
pact characteristic data was studied in Ref. 6], where 
the profile of the characteristic scalar field was logarith- 
mic in advanced time. These data preserve the picture 
of the singularity inside spherical charged black hole: the 
singularity is a union of spacelike and null sectors. The 
causal structure of spacetime is again described by the 
lower diagram of Fig.^ In this paper we show that other 
families of characteristic data may lead to dramatically 
different causal structures. 

We write the spherically-symmetric metric in double- 
null coordinates in the form 

ds^ = _2e2'^("''') du dv + (u, v) d9? (1) 

where dfi^ is the line element on the unit two-sphere. As 
the source term for the Einstein equations, we take the 
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FIG. 2: Spacetime diagrams with variable A. Lower diagram: 
Spacetime diagram for a > Qcrit. Upper diagram: Spacetime 
diagram for a < Ocrit . In both cases spacelike singularities are 
depicted by a bold solid line. The outgoing null ray (dashed 
line) at m, intersects with the spacelike singularity at a corre- 
sponding afiine parameter (along uo) 7* (dotted lines). With 
lower value of A,the event horizon (EH) moves inward, and 
the spacelike singularity occurs later. 

contributions of both the scalar field $ and the (source- 
less) spherical electric field (see for details). The dy- 
namical equations are the scalar field equation □$ = 
and the Einstein equations, which reduce to 

+ - {r,u'^,v + r,^<^,u) = (2) 
r 

r,uv + '- + ^ 1 - — = (3 

r 2r \ J 

and 

'-.^ - ^ - ^ (l - 2^ j + ^.^^.^ = ■ (4) 

These equations are supplemented by the two constraint 
equations 

r„„ - 2cr,„r„ -hr($,„)2 = (5) 

-2CT^r^, -f r($,^,)2 = 0. (6) 

The initial value problem for these equations and our 
choice of gauge for the coordinates are described in 
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FIG. 3: AfRne parameter at uo 7* at which a spacehke singu- 
larity is intersected along an outgoing null ray m* as a function 
of A. Here, a = 0.3. Numerical data are depicted by circles, 
and the solid line is proportional to j4~* °°. 



FIG. 4: Afhne parameter at uq 7, at which a spacelike singu- 
larity is intersected along an outgoing null ray u, as a function 
of A. Here, a = 0.1. Numerical data are depicted by circles. 



Ref. la]. We choose a gauge in which r = w on the out- 
going segment of the characteristic hypersurface u = uq 
and r = To + uTuo on the ingoing segment, v = vq. We 
then take the characteristic initial data to have a profile 
of <I> = ^(w" — Vq) on u = Uo, and ^ = on v — vq. 
The characteristic value problem is then satisfied by tak- 
ing a = aA^iv^"' - vl°^)l\ - (1/2) In 2 on m = uq and 
(T = — (l/2)ln2 on V — vq. The parameter r„o is deter- 
mined by the initial mass of the black hole Minitiai, its 
fixed charge Q, and r(uQ,uo). Here, u,v are (twice the 
standard) retarded and advanced times, respectively. Be- 
fore the onset of the perturbations {v < vq) the geometry 
is that of a RN black hole. We then vary the exponent 
a, and study the occurrence of a null singularity with 
different values of a. 

For a given value of a, the most straightforward way 
to establish the occurrence of a null portion of the sin- 
gularity is to consider a given value for A, and then 
probe the fields along various outgoing null rays. In- 
deed, that was the method used for finding the null sin- 
gularity, and the transition from a null singularity to a 
spacelike singularity (at later retarded times) in Ref. . 
In contrast, here we wish to establish the inexistence of 
null singularities for classes of families of characteristic 
data. As numerically it is not easy to locate outgoing 
null rays at arbitrary values of retarded time, and as 
the level of required fine tuning is unbounded, we un- 
dertake here a different approach: We fix the retarded- 
time value of an outgoing null ray (using the geometry 
before the onset of the perturbations, i.e., for v < vq) 
at u = > Mo, and vary the amplitude A (for fixed 
values of a). (We consider here positive values of A.) 
That way, we consider a family of spacetimes (sharing 



a common a). In order to compare advanced times be- 
tween different spacetimes, we reexpress advanced time 
in terms of the affine parameter 7 along u = uq. It can 
be readily shown that along uq (and for fixed angular 
coordinates), 7(1;) = vM[l/{2a), 1 + l/{2a),aA'^v'^"' /2], 
where M{a, b, z) is Kummer's function. Along u*, 7(f) is 
of course not an affinc parameter, but wc can express v 
along u■^, in terms of 7 using ingoing radial null geodesies. 
We denote by 7* the value of 7 at which a spacelike sin- 
gularity is intersected hy u — u^. We remark, that these 
results show that u = uo is geodesically complete. Notice 
that 7(w) is a monotonically increasing function. (We do 
not have problems using u when comparing between dif- 
ferent spacetimes, because u is defined in terms of the 
RN geometry that preceeds the onset of the perturba- 
tions.) With smaller amplitude A, the intersection of the 
null ray with the singularity occurs at a later value of 7* . 
(See Fig. 13) That is, 7, is a decreasing function of A. If 
a null singularity exists, there is a finite, nonzero value 

transition, SUCh that 7* ^ CX3 aS A ^ ^iansition 7^ 0. If 

no portion of the singularity is null, 7* — > 00 as ^ — ^ 0+ . 
[The case A — corresponds to the exact RN solution 
(no perturbations).] That is, the singularity is spacelike 
for arbitrarily weak perturbations, and the causal struc- 
ture of spacetime is described by the upper diagram in 
Fig. in 

The results described below refer to the following 
choice of parameters: Q = 0.95, Afinitiai = 1, r(Mo,wo) = 
5, and = 21.4. Qualitatively simlilar results were ob- 
tained for other choices of the parameters. The numerical 
code is based on the code in Ref. ^ . We tested the code 
and found that it is stable and second-order convergent. 

Next, we set a = 0.3. Our results are shown in Fig.O 
which displays 7* vs. A. For small values of A, 7* (A) cx 
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FIG. 5: The exponent /3 as a function of a. The circles are 
the data points, and the solid curve is proportional to a~^'^^ . 



, where jS = — 8.00±0.01. Extrapolating this behavior 
to smaller values of A we conclude that 7* (A) ^ cxj as 
A ^ 0^. That is, for this choice of a wc find that only 
a spacelike singularity exists, and that there is no null 
sector to the singularity. In Ref. it was shown that 
a spacelike singularity can exist inside spherical charged 
black holes. Here, we show that inside some spherical 
charged black holes the singularity is only spacelike. 

Figure 21 shows analogous results with a = 0.1. Here, 
we find (using Richardson's deferred approach to the 
limit) that 7*(^) ^ 00 as A ^ Atansition = 4.26 ±0.01. 
That is, for A > ^transition the outgoing null ray inter- 
sects with a spacelike singularity, and for A < ^transition 
it intersects with a null singularity. 

In between these two values of a we expect that a 
critical value exists, which marks the transition between 
black holes whose singularity is composite of two types 
(spacelike and null) , and black holes whose singularity is 
an all-encompassing spacelike singularity. 

We return now to the supercritical case in which only a 
spacelike singularity exists. The exponent /3 is a function 
of a. Figure shows /3(a). The exponent /3 is a mono- 
tonically increasing function of a, and approaching (from 
above) the critical value acrit its slope becomes steeper. 
Fitting this curve, we find that f3{a) oc q,-1-3i±0.07^ 
tice, that (3 is only defined for a > acrit- 

We find that the critical exponent is acrit — 0.12±0.02. 
The corresponding value of /3 is /3crit = —26.9 ± 3.0. 
That is, for values of a smaller than acrit we find that 
both spacelike and null sectors coexist. For such values 



the resulting spacetime has the same causal structure as 
spacetimes which evolve from perturbations with com- 
pact support. In particular, the survival of a null sec- 
tor leaves open the possibility of harmless traversal of 
the null singularity. In contrast, for values of a larger 
than acrit no null sector survives. The spacelike singu- 
larity is all-encompassing, deformationally strong, and 
completely blocks the "tunnel" inside the black hole. 

The numerical value we find for acrit is not expected 
to be universal, but instead be model dependent. In par- 
ticular, we expect it to depend on the ratio Q / Minitm\, 
and the type of field used. In particular, we expect acrit 
to increase with Q/Minitiai- Our main result, namely 
that there are families of characteristic data for which 
the singularity inside a spherical charged black hole is 
an all-encompassing spacelike singularity, is expected to 
be held in general. It is an open question whether our 
results hold also for spinning black holes. The question 
of whether the transition we find resembles mathemat- 
ically the critical phase transitions of statistical physics 
(like the threshold of black hole formation 9] ) awaits fur- 
ther investigation. If so, a possible candidate for a scaling 
relation could perhaps be the affine length of the null sec- 
tor of the singularity as a function of acrit — a. Finally, 
we comment that realistic black holes are always irradi- 
ated by photons which originate from the relic cosmic 
background radiation. This paper then is relevant to the 
question of whether these photons completely block the 
"tunnel" inside astrophysical black holes. The answer 
to this question certainly depends on the cosmological 
model. 

I thank Chris Beetle and Richard Price for discussions. 
This research was supported by the National Science 
Foundation through Grant No. PHY-9734871. 
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